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Nonuniform lossless or lossy transmission lines or layers used as broadband matching 
or absorbing devices are studied. 

When the refraction index, n(x), and the characteristic impedance Z (x), are given, the 
reflection spectrum, p (??) =po(4x/X), for x = can always be computed by solving numerically 
a Riccati differential equation (RDE). (Analysis) 

Conversely, not only for n = const [Bolinder, 1950, 1956] but also for n(x) real and m = Mo, 
a tapered transformer can be synthesized starting from a given pq(v) spectrum by using 
Fourier transform techniques. (Synthesis) 

For broadband absorbers, the synthesization procedure can be approximately applied, 
under certain conditions, to only the part of the spectrum which represents the reflection of 
the matched (lossy) line. 

1. Introduction 

The first part of this paper deals with the problem of matching two uniform transmission 
lines (or layers) by means of a nonuniform line (or layer) in such a way as to "minimize the 
reflections." 

A variational solution of the problem does not seem to exist, and, therefore, instead of 
trying to minimize the average value of the reflection over a given frequency band, we will 
solve the problem of keeping the reflection lower than a given value over the largest possible 
frequency band, for a given length of the matching line. 

To find an acceptable solution for this problem, two methods can be followed. The 
first one (Analysis) consists in varying tentatively the parameters of the line in such a way 
as to obtain a satisfactory reflection spectrum, i.e., reflection coefficient versus wavelength. 
According to the second method (Synthesis), we start from a favorable reflection spectrum and 
find the corresponding law of variation of the parameters of the line. 

The second method is made possible by the fact that, in a number of cases, the line 
parameters can be obtained operating on the inverse Fourier Transform of the reflection 
spectrum. 

In the last paragraphs of this paper (par. 6 and 7) some tentative techniques are described 
to utilize even in the case of absorption the synthesis method, which is rigorously appliablc only 
in the case of matching. 

Let us now recall the basic equations to which recourse will be made. 

For monochromatic waves, the two wave equations describing the propagation of em 
waves along an inhomogenous line or a stratified medium can be condensed into the single 
Riccati differential equation (RDE) 

Z'+yZ -y^=0 (1.1) 



1 The research reported in this document has been sponsored by the European Office, Office of Aerospace Research, USAF, under Contracts 
A F61 (052) -36 and 589. 
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where 

Z— input impedance at abscissa x, 

Z' = its derivative with respect to x, 

Z = intrinsic impedance of the line or layer at abscissa x, 

2tt 
y =j— n (x) — propagation constant, and 

A 

n (x) = Tii (#) —j^2 (x) = relative refraction index. 

Equation (1.1) can be written in terms of the reflection coefficient p(#), which can be de- 
fined in the "Schelkunoff form" 

Z(x)—Z (x) . 

pM =Z(x)+Z (xy (L2) 



obtaining the following RDE : 



p , -2tp+^(1-p 2 )=0. (1.3) 



Let us note here that (1.1) and (1.3) are valid for layered media not only under normal 
incidence, but also under ah conditions of incidence and polarization, provided the normal 
equivalent wave is considered. 

If we multiply both p. r (x) and e r (x) by a constant C, the first and last terms in (1.3) remain 

4-7J- 

unchanged, while the second becomes —j~—Cn(x)p(x). It follows that we find the same 

A 

reflection coefficient p(z), for the wavelength X, that we found previously for the wavelength 
\jC; that is, on the average, the reflection coefficient is lowered. This property cannot be, in 
general, easily utilized in the design of inhomogeneous layers, since high values of p. T are usually 
accompanied by still higher values of e r . The above remark shows, however, that an increase, 
even if gradual, of p r (winch is generally a complex quantity), increases the effective length of 
the line, and lowers, consequently, p(x). 

When ji=const, and only e varies with x, a single function is sufficient to describe the 
inhomogeneous line: either y(x), or Z (x), or n(x). Equation (1.3) can then be written in the 
form 

p'-:hnp-^(l-p 2 ) = 0, (1.4) 

where 

77 = 47r/A. 

2. Analysis of an Em Line or Layer 

The problem of analyzing the behavior of an em line or layer can be summarized as follows: 
For given y(x) or n(x) and Z (x) functions (i.e., for given e,(x) = e rl (x)— je r2 (x) and ii r (x) = iJL rl (x) — 
jn 7 2(x) functions) and length L of the layer, find the reflection coefficient at abscissa x=0, 
p(0, v) = Pq(v), as a solution of (1.3) or (1.1) plus (1.2), with a given limit condition p(rj, L). 
Obviously, when pt=const it is sufficient to start from only one of the previous functions, e.g., 
from the refraction index n(x). 

In the case of matching, p(r), J,) = 0; in some cases of absorption, p(rj, L) = — 1, i.e., the line 
is short circuited at the far end. It must be pointed out that (1.1) is more general than (1.3), 
since it allows an initial discontinuity. If there is such a discontinuity, only if both this one 
and po are very small, the overall reflection coefficient can be approximately calculated by 
adding the two reflections. 

Generally, it is not possible to have analytic solutions of (1.1) and /or (1.3), and it is neces- 
sary to carry out computations by means of digital or analog computers. 
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Equation (1.4) can bo linearized as follows: 



p'— jynp-7fo =0, 



(2.1) 



provided that |p(x)| 2 <<l everywhere. This condition is generally fulfilled in a matching 
device, at least in the rj region of low reflection. 

In the case of matching and for p r =l, ten different real n(x) functions were tested, and the 
corresponding p Q (rj) computed from (2.1), by means of a PACE-TR-10 analog computer. 

Assuming n(0)=l, the "transformation rates" n(L)/n(0) = 2, 5, and 10, were considered. 
It was shown [Franceschetti, 1962] that, for n(0)^l, the results are still valid, provided that 
we refer to wavelengths measured in the medium of refraction index n(0), i.e., to the wavelength 
\/n(0). The scattering matrix (S) of the lossless matching junction can be easily deduced 
starting from the knowledge of p (rj) = S n , computed by means of the linearized (2.1). Actually 
it is possible to express the coefficient S 22 in terms of p (??), and, by means of the matrix equation 
(S)(S) *=(/), to calculate (S), obtaining the following expression [Franceschetti, 1902]: 



(S) = 



pM Vl— |Po(*?)| 2 -exp — j~ I n{x)dx 

Vl — |p (*?)| 2 -exp I — j I I n(x)dx -pJW-exp —jrj n(x)dx 



(2.2) 



where p$(ri) is the complex conjugate of p (??). 

The results of the p Q (rj) calculations, in the ten different cases considered here, are referred 
to under table 1, where the ratio £=A/Z up to which |po(r?)| % is less than 10 percent is given 
for the ten n(x) functions, and the three transformation rates 2, 5, and 10. This means that, 
for wavelengths A<££, the matching junctions under consideration give a power reflection 
less than 1 percent. It is apparent from this table that, assuming |p (^)|< 10 percent, a 
generally "optimizing" n(x) function cannot be defined, since the behavior of each function 
depends on the transformation ratio. For example, the best matching function seems to be 
7i(z)=exp [kx/L], for n(L)=2 (£=3.85); and ^(x)=exp [&(z/Z) 3 ], for n(L)=5 (£=2.82), and 
n(L) = 10 (£=2.18). 



Table 1. Values o/£ = \/L such that for X<L^,|p (L/X)| < 10 percent (matching case) 





exp 

[kx[L] 


exp 


exp 
[k(x/L)»] 


1+kx/L 


l+k(x/L)* 


\+k(xlL)* 


1 


exp 
[ekx/L-i] 


cosh 

[kx/L] 


1+k 


Functions n(x) 


1-A- x/L 


(1-COSttX/L) 


n(L)=2 
n(L)=5 

?7(X)=10 


3.85 
1.36 
1.44 


2.60 
2.0 
1.89 


2.04 
2.82 
2.18 


3.78 
<1.0 
<1.0 


3.35 

2.50 
1.61 


2.62 
1.99 
1.46 


3.57 
1.58 
1.45 


3.64 
2.41 
1.80 


3.08 
2.41 

1.77 


2.94 
1.33 
<1.0 



It must be pointed out that values of n up to 8 are easily attainable b\^ means of metal 
powders suspended in paraffin wax [Kelly et al., 1953]. 

Tn the case of absorption, ten different nonmagnetic and magnetic functions n(x) were 
tested, and their behavior computed, solving (1.1) and (1.2) (i.e., the nonlinearized equation) 
by means of the IBM 1620 digital computer of the Faculty of Science of the University of 
Naples. An initial discontinuity was often assumed. 

The results are summarized under tables 2 and 3, where the power reflection coefficient 
|p | 2 % is given for several values of the parameter X/L. 
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Table 2. Power reflection coefficient |p (L/X)| 2 percent, for several values o/X/L (absorbing case) 
e=Aexp[2(m-j<p)x/L] n=Bex\) [2(n-jt)xJL] 



A.m. <p 
B.nA 


1.0 1.0 0.7 
1.0 0.0 0.0 


1.5 1.0 0.7 
1.0 0.0 0.0 


1.5 1.5 0.6 
1.0 0.6 0.6 


1.5 1.0 0.7 
1.0 1.0 0.6 


1.0 1.0 0.7 
1.0 0.5 0.5 


X/L = 1.0 
2.0 
3.0 
4 
5.0 
6.0 
7.0 


0.30 
1.20 
2.30 

.80 

8.2 

23.8 


0.75 
1.0 
2.70 
2.90 
.50 
3.75 
14.8 
30.16 






0. 054 


2.31 
2.92 
3.07 
4.12 
4.95 
5.60 
6.04 


0.61 
.60 

.82 
.80 
.80 
.72 
.59 


.29 
.45 
.59 
.66 
1.06 
1.99 


8.0 




3.45 









Table 3. Power reflection coefficient |p (L/A)| 2 percent, for several values of X/L (absorbing case) 

x X 
e=A [cos 2<p - -j — sin 2<pXJL]- cxp [2m x/L] /x=Bexp [2{n-j\p) x/L] 



A. m.<p 


1.5 1.0 0.7 


1.5 1.5 0.3 


1.5 1.5 0.3 


1.5 1.5 0.3 


1.5 1.5 0.3 


B. m.t 


1.0 0.0 0.0 


1.0 0.4 0.4 


1.0 0.4 0.4 


1.0 0.6 0.4 


1.0 0.6 0.4 


T 


5 


3 


5 


3 


5 


X/Z = l 


1.70 


1.83 


1.56 


1.75 


1. 58 


2 


.46 


2.13 


2.24 


1.91 


2.01 


3 


5.90 


2.99 


3.26 


2.92 


2.12 


4 


2.90 


3.87 


3.76 


3.12 


3.34 


5 


.50 


5.46 


5.51 


3.94 


3.76 


6 


4.50 


7.47 


7.52 


5.36 


5.20 


7 


12.4 


9.35 


8.21 


7.03 


6.45 


8 


22.3 


11.0 


8.26 


8.74 


7.03 



The values of the parameters are listed in the tables under each case. 

In table 3 the conditions for physical realizability (i.e., the Kramers-Kronig equations) 
[Latmiral et al., 1961] have been taken into account for the electric losses (purely conductive 
electric losses have been assumed) . 

Cases 1, 2, and 6 are nonmagnetic; besides, for x— >L, they present values of the rate 
e 2 /€i so high as to be difficult to obtain practically. 

Cases 3, 4, and 5, and particularly the last, approach too closely to Heaviside conditions; 
they assume values of \x r practically unattainable in the microwave range. 

Cases 7, 8, 9, and 10 are probably those which may be realized without excessive diffi- 
culties; the electric losses are purely ohmic, and the rate <r/coei seems to be everywhere realizable. 

3. Problem of the Synthesis 

The problem of synthesizing an em line (or layer) can be summarized as follows: 

Starting from a given (generally complex) reflection spectrum p (r)) at the beginning of 
the line of length L and from given limit conditions, find the two corresponding y(x) (or n(x)) 
and Zq{x) functions. If /^^const, only one real function, e.g., n(x)> must be found. In the 
case of matching, obviously, the limit conditions are the values of n(0) and n(L). 

In the following sections it will be shown that the synthesizing problem for the match- 
ing case, at least for nonmagnetic lines, can be completely solved by means of the Fourier 
transform techniques. 

On the contrary, in the case of absorption (except for some special cases), only the follow- 
ing problem can be solved: for a given 2 (^) function, and length L of the line, find a com- 
plex function n(x) such that its reflection coefficient p (rj) is more favorable, on the average, 
than the given 2 (*?). 

4. Synthesis in the Matching Case 
As was pointed out in section 2, in the case of matching the linearized RDE, 



P '-2 7P +^=0 



2Z n 



(4.1) 



may be considered instead of (1.3). 



334 



Solving (4.1) for p(rj, L)=Q, we have 

PoM= I 2Z exp | ~ /r? I ^0*)^ P x - ( 4 - 2 ) 

When n= const (this is the case of TEM guided waves), and, e.g., n=l, (4.2) becomes 

c L z r 

poM= I 2Z exp [~^ rfx - ( 4 - 3 ) 

Equation (4.3) shows that p Q (rj) is the Fourier transform [Bolinder, 1950, 1956] of the 

Z' 

function (of x) -^d which is zero out of the interval (0, L) in which the inhomogeneity is con- 

fined, and which in turn can be represented as an inverse Fourier transform as follows: 

z' i r +o ° 

2Z=2~A _ PoM explfax] dv- (4.4) 



Z' 

The fact that p (rj) and ^Mr are a transform pair provides the argument for making the 

synthesis of a broadband matching line with 7&=const. 

Z' 

For example, if we choose for p (rj) the well known "sampling function," we must have -%- 

I/jo 

constant in the interval (0, L) (and zero out of it). Analogously, if p (rj) is the square of the 
sampling function (which is a very favorable "spectrum" of reflection), ■=£- will become a 

"triangular pulse" in the interval (0, L) [Bolinder, 1950, 19561. 

When ;u=const, we can get a transform pair even when n is not a constant, but a real 
function of x. 

The constancy of p. allows us to start with (2.1) instead of (4.1), so that (4.2) becomes 

Po(i?) = — j ^expl — jy) j n(x)dx \ dx. (4.5) 

As n(x) is always positive, the integral function 

y=y(x)= n{x)dx (4.6) 

is always increasing and admits, consequently, the inverse function x=x(y). 
Let us introduce in (4.5) the new variable of integration y. We obtain 



PoM 



-~fS§H «*-'-*"*• <«> 



where y L = y(L), and the derivative is with respect to x. 
As (4.7) shows, po(ry) and 

, x n'[x(y)] ( . 

v^-mxJyT] (4 - 8) 

is a Fourier transform pair. 

As a matter of fact, even when y L is finite, the function p(y) is zero out of the (0, y L ) in- 
terval and therefore the integration can always be understood as extending along the whole 
positive y axis. 

If the po(ri) function is given, the <p(y) function result determined by means of the inverse 
transform and the n(x), which is the datum of practical interest, can be obtained by solving 
the integral DE 

n'(s) 



2n 2 ( 

672865—63 6 335 



^=p(J\(z)ix), (4.9) 



which becomes a DE of the second order in the unknown y= \ n(x)dx and, precisely, 

J o 



-f^=^)- 



(4.10) 



Let us start from a family of spectra p (r], A, y L ), whose inverse transforms are zero out of a 
finite interval; the parameter y L represents the extent of said interval. 
Equation (4.10) becomes 



V 

-fyT2=<p(y,A,y L ), 



(4.11) 



where the derivatives are with respect to x. 

This is a DE of the second order containing two parameters (A, y L ). 

Integrating the above DE with the two initial conditions t/(0)=0 and y' (0)=?i(0), we 
get a solution which contains the parameters A } y L : 



y=y(x,A,y L ). 
The parameters must be chosen in such a way that 

y L =y{L,A,y L ), 
n(L)=y'(L, A,y L ). 



(4.12) 

(4.13) 
(4.14) 



Obviously, for ij— >0, i.e., for very long wavelengths, the reflection coefficient p (r]) must 
approach the value [^(0)— n(L)\ ■ [n(0)+n(Z)] _1 . However, it [must be pointed out that we 
do not deal here with the complete RDE, but with the linearized one, and that the latter does 
not require the above condition, which will be satisfied only for low "transformation rates." 




Figure 1. n(x) functions referring to table 5, 1st 
row. 

The various curves are related to different values of the transforma- 
tion ratio n(L)/n(0) = n(L)/l=n(L). 
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In some special cases, the problem can be solved analytically [Latmiral et al., 1962]. 
For example, when p (r?) is of the type of the "sampling function," the corresponding n(x) 
function becomes of the type proposed by Jacobs [1958] on the basis of empirical arguments. 

Analogously, when p (rj) is an exponential integral (which is a very convenient type of 
reflection spectrum), the corresponding n(x) function becomes the well-known exponential 
one [Latmiral et al., 1962]. 

In general, however, the problem admits of only numerical solutions, and it is necessary 
to solve (4.11) under conditions (4.13)-(4.14) by means of analog or digital computers. 

Five different types of <p(y) functions (to which favorable spectra of reflection are related) 
were chosen, and the corresponding n(x) functions computed, for several rates of transformation, 
by means of the above procedure. 

The n(x) functions are traced in figures 1 through 5, under each of which the corresponding 
<p(y) function is referred. The first four cases were computed by means of an IBM 1620 digital 
computer, and the last one by means of a PACE-TR-10 analog computer. 

Besides, in table 4, the ratio f=X/Z, up to which the reflection coefficient |p (*?)| is less 
than 10 percent, is given for the computed n(x) functions and three transformation rates 
(2, 5, and 10). The complete p (^) function can be easily computed as the Fourier transform 
of the corresponding <p(y) functions, and are referred, together with the numerical values of 
the parameters, in table 5. 

In all computations we have chosen n(0) = l, but this assumption does not cause any loss 
of generality (see sec. 2). 

As in the case of analysis (see sec. 2), it is apparent that a generally "optimizing" function 
n(x) cannot be defined. It is interesting to point out the very good results obtained for high 
"transformation rates." For example, for n(L)/n(0) = 10, the n(x) function synthesized starting 
from the function <p(y)=A sin [iry/y L ], presents a value of {=3.80, i.e., wavelengths up to 
about four times the length L of the matching junction are almost completely transmitted. 




n(x) 




Figure 2. n(x) fund ions referring to table 5, 2d 
row. 



Figure 3. n(x) functions referring to table 5, 
row. 
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Figure 4. n(x) functions referring to table 5, 4th Figure 5. n(x) functions referring to table 5, 5th 

row. row. 

Table 4. Values of £ = X/L, such that for \^L £, | po(Lj\) | 2 — 10 percent (matching case) 



Functions <p (y) 


A[yyL-y 2 ] 


A[yy L -y 2 ] 2 


.4(1 -cos 2ir ylm) 


A sin 7T yfyL 


Triangular pulse 
between 0-?//. 


n (L) = 2 
n (L)=5 
v (L) = 10 


2.90 
3.20 
3.70 


2.08 
2.£0 
2.80 


2.02 
2.50 
2.63 


2.90 
3.20 
3.80 


2.27 

2.63 
3.03 



Table 5. Matching case: Functions <p(y), P0O7), and numerical values of the parameters 



Fig. 



<p(y) 

V<\I<IIL 



Pq(v) 



n(L) 
(exact 
values) 



A(yy L -y*) 



A{yyi-yW 



l(l -eos|-0 



A sin - 



Triangular pulse of 
height A 



2/1 T2 . 77.VL y]VL~\ f • Vl~\ 



AyL 



J b sm a T cos V ""' sm "2" J • cxp L - Jv 2 J 



2 



IGz*]*} [£-}¥] I J 



•M] 



LkHjJ Lfc~*j7j J 



P[^]'-[> ? l] 



Ai/l fsin riyLn 

8~ 



2.05 
2.94 
5.05 
6.99 
11.0 

2.01 
2.97 
4.95 
7.20 
9.85 

2.01 
2.99 
4.97 
7.01 
10.07 

1.99 
3.10 
5.20 
6.68 
10.90 

2.0 
3.0 
5.0 
7.0 
10.0 



-.61 
-.55 
-.46 

-1.95 
-1.56 
-0.75 
-.72 
-.53 

-.25 
-.32 
-.40 
-.49 
-.52 

-.36 
-.54 
-.65 
-.68 

-.75 

-.50 
-.704 



1.5 
1.7 
2.0 
2.2 
2.5 

1.4 
1.6 

2.0 
2.1 
2.3 

1.4 
1.7 
2.0 
2.0 
2.20 

1.50 
1.68 
2.00 
2.20 
2.50 

1.35 
1.60 
1.86 
2.11 

2.28 
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5. Quarter Wavelength Law 

Let us now define the interval A-q according to the following relation (analogous to the one 
defining the "radius of gyration" in mechanics, or the root-mean-square error in error theory) : 



where 

/ 



(Av) 2 = 7 r\' 2 \poM\ 2 dv, (5.1) 

* ./ — 00 

r° \PoW\ 2 dr,. (5.2) 

7 — CO 



Let Ay be the analogous interval for the function <p(y). 

It is apparent that the two functions p (i?) and <p(y) will be noticeably different from 
zero only inside the two intervals A77 and Ay. 

The product of the intervals Arj and Ay cannot be less than a given positive value K, which, 
according to the adopted criterion of the square mean, equals 2x [Persico, 1950]: 

Ar).Ay>2ir. (5.3) 

If we refer to the positive half interval, A^=47r/X max , where X max is the longest wavelength 
almost completely transmitted, and if we put Ay^y L =7iL, where n is the mean value of n(x) 
in the interval (0, L), (5.3) becomes 

2^-nL>2v X max <4rLL. (5.4) 

A m ax 

When n=const=l, the two Fourier related functions are p (r]) and Z /2Z (see sec. 4), 
and inequality (5.4) becomes 

X max <4Z. (5.5) 

This last inequality shows that, for a given length L (in air), the matching line becomes 
ineffective for wavelengths beyond 4i, whatever the function Z' Q /2Z may be. 

On the contrary, when n is a function of x, at least in nonmagnetic cases, an improvement 
can be expected according to (5.4). 

For further improvements due to magnetic properties, see section 1. 

6. An Approximate Solution of the RDE in the Case of Absorption 

For the synthesis of a broadband termination, an empirical method may consist of adding 
an arbitrary imaginary part to a purely real n(x) function, found according to the above tech- 
niques (see sec. 4), in such a way as to make the phase vary linearly with x; a numerical checking 
is then obviously necessary. Another way may consist of planning an absorbing device com- 
posed by a broadband transformer with a high n(L)fn(0) rate (table 4, n(L)/n(0) = 10), followed 
by a thin "Heaviside" absorber, probably realizable without excessive difficulties in these 
conditions. 

We will now find an approximate solution of the RDE in the case of absorption, with the 
limit condition p(t?, L) = — 1. 

As is well known, a RDE can be transformed into a Bernouilli one, once a particular 
integral is known. We may consider as such the solution of the RDE for the condition 
which corresponds to perfect matching at the far end, i.e., to p(rj, L)=0. 

Calling <t(x) such a solution, the substitution r=p—a transforms (1.3) into a Bernouilli DE 
in which the further substitution y=r~ l yields to a linear equation whose general integral is the 
following : 

,=exp [£(*+. f ) &] { ~\ £ I exp [-f (*+. § *}*+<7} • (0.1) 
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As y=(p—a) l must become ( — 1—0) 1 = — 1 at the end of the shorted line, C must equal 
— 1 . Consequently, we can write : 

ex P r-rv2 T +«rf)^] 

1 -2j.^; esp L"J. v y+a zj (h \' h 

For x=0 we have 

expf-pfVl-fffVzl 

/ \ /x L Jo \ A)/ J / P Q \ 

PoU) = *oU) 1 /» L Z / F /. L / — y7\ n («.«$) 

'-»/. I esp [-X 0+'i>*> 

As o-(z) vanishes at the end of the line, contrary to p(x) which equals —1 at the same point, 
the RDE for <r(x) can be linearized, neglecting |cr(x) | 2 in comparison with unity; consequently, 
we can consider a(x) as the solution of the following linear DE : 

Multiplying (6.4) by a and neglecting a 2 , we have 

do- 2 =-§ <rdx. (6.5) 

By integration we obtain 

- F ^ dx=-o>(x). (6.6) 

Consequently 

exp |~- p a ^ dafUexp [-o*(x)]^l-o*(z), (6.7) 

and, analogously 

exp [_— J o <r^(faj=exp[-a5(ij)]^l-aj(i|). (6.8) 

Substituting (6.7) and (6.8) into (6.3) we obtain 

[1-ogfo)] exp [-2 J%fe] 

Po(^7) = ^o(^) 1 /»& Z / T^ TTl =j (6.9) 

1— JJ o g-[l-c7 2 W]-exp[-2J^ T ^J.r/.r 

Neglecting c 2 we get the following approximate value for p (rj) : 

exp —2 7(fo; exp —2 7<fa 

Po(i?) = <ro(*?) , /»& 7/ f^"^ — n — =(r o(^)— - 1 \-\\ ' (6.10) 

where 

^(n)=+g jjf exp t" 2 /^] ■ ^ (6 - n) 

and <t l (L) is the reflection coefficient in x=L for a wave which proceeds from right to left in the 
layer matched at x=L [Latmiral et al., 1962]. 

When |a : L (r7)|<<l, (6.10) can be further simplified as follows: 

PoW^oW-exp \j-2J o ydxj, (6.12) 
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even when 7 is complex, and n is not a constant but a complex function of x. 

It must be pointed out that (6.12) is the solution of the RDE (1.3) linearized by neglecting 
|p 2 | against unity and solved with the condition p(L) = — 1. 

This simplification is, however, possible only as a consequence of the considerations which 
have led to (6.10) by splitting the solution of the RDE into two parts by means of the particular 
integral a(x), which vanishes for x=L. 

Otherwise, the two assumptions |p| 2 <<l and p(L) = — 1 would have been almost un- 
justifiable. 

7. Special Techniques for the Synthesis in the Case of Absorption 

The splitting of the solution of the RDE into two parts (6.10), the first of which refers to 
the condition a(L)=0, allows us to extend the transform techniques (under the conditions 
pointed out in sec. 6) even to the case where the layer is not matched at the far end and par- 
ticularly to the case where p(Z) = — 1. Equation (6.10) shows, as a matter of fact, that to 

minimize \po(rj)\ is equivalent to minimizing |o- (t;)|, provided that i(x)=] ~ n(x) is chosen in a 
class of functions which satisfies the condition 

(Real part of f%k)==| ( L fh(z)dx>H, (7.1) 

where H is a given positive quantity. 

According to (6.12), this means that the wave reflected from the metal must not exceed a 
reasonably low value. 

Obviously, in the case p(L) = — 1, n x {L) andr& 2 (£) are not given, and only fti(0), n 2 (0) and 

n 2 dx (see 7.1) are given. 

The above considerations allow us to extend, at least in some special cases, the synthesizing 
procedure referred to under section 4, to the a (rj) function [Latmiral et al., 1962]. 

For example, when the refraction index is of the type n(x)=rii—jK/ri } with rii and K real 
positive constants with respect to both 77 and x, the Fourier transform pair (see sec. 4) becomes 
Po(tj) and (Zq exp [— Kx])/2Z . This case occurs, e.g., when very low ohmic losses are present 
(o-<^<^coe), and e and a do not depend on x. 

Another case in which the extension is possible is n = a l —jn 2 = const, with n x and n 2 both 
independent of x and 77. p (ri) are then the values that the Laplace transform of the function 
Z'J2Z takes on the straight line of complex equation 

p=(n 2 +jni)r} } (7.2) 

with 77 in the interval (0, 00). 

When the refraction index is x-dependent, the synthesizing procedure is still applicable 
in some special cases. 

For example, when the refraction index is of the type n(x)=ni(x)—jK/ri, with inconstant 
with respect to both 77 and x, (4.11) becomes, under the hypothesis that e rmax is not too high 
(for example, € rmax <16), 

-f^e- K *=a (y)- (7.3) 

This DE must be solved according to the method explained in secti on 4 . The above 
case occurs when very low ohmic losses are present, and the quantity <r(x)/^e r (x) (proportional 
to the loss angle) is constant with respect to x. Another case is the following: 

n(x) = (n l -jn 2 )f(x), (7.4) 

with rii and n 2 real constants with respect to both 77 and x. 
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Let 



V= f(x)dx; y L = f(x)dx; Z [x(y)] = Z (y). (7.5) 

Jo Jo 



The values of p (rj) are those which the Laplace transform of the function Zq/2Z takes 
on the straight line of complex equation [Latmiral et al., 1962] 

P=(ri2+jn 1 )r 1 . (7.6) 

8. General Techniques for the Synthesis in the Case of Absorption 

Let us now consider the general problem of synthesizing a nonmagnetic absorber; we will 
solve this problem by means of an approximate method. 

It must be pointed out that, in the absorption problem, not four but only three conditions, 



wi(0) 



, 7i 2 (0), and n 2 dx, are given in a nonmagnetic absorber. The following approximate 

Jo 

procedure, based on transform techniques, may be helpful. For /ti= const and small values of 
p(x), the linearized (1.3) (see sec. 6) gives: 

<ToO?)=~ 2 J nl-jn 2 eXp L~ V Jo ^ J ' 6XP L^^Jo ^J **' ^'^ 

Let us compare the above equation with the following : 

s.««4 f ^t; exp [-*M dx - (8 - 2) 

If some exceptional cases are excluded, the value of 2 (*?) are, on the average, less than the 
values of a (r)). 

Consequently, if rii(x) and n 2 (x) are such as to keep S small in a given interval of rj, the 
same functions can give a convenient <r . 

Introducing the new variable of integration 

V=fni(3)<fe, (8.3) 

So(*)=— g J o V [In (n 1 -jn 2 )]> exp [—jny]dy. (8.4) 

Thus 2 (t?) is a Fourier transform. 

Once 2 (r?) has been chosen, and indicating its inverse transform by (p(y)= z <pi(y)—j(P2(y), 
we have 

2dy ^ ln C^i — i^2)]=^i(y) — jV2(y), (8.5) 

| (In [%(2/)-> 2 (7/)]-ln [MO)-jn a (0)])=£^(y)dy-jj^(y)dy; (8.6) 

In [n, (y)-jn 2 (y) =( 2 /Vi (y)dy+h^n>M+n 2 2 (0)) 

-j(2jV 2 (y)^+arctg^)=^(2/)^j^ 2 (y). (8.7) 
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i.e. 



As 7h(0) and n 2 (0) are given, the values of n x (y) and n 2 (y) can be obtained in terms of y, 
solving (8.7), i.e., 



^/nl(y)+n 2 2 (y)-=cxp [%(y)]; arctg 



My) _ 
My)' 



^2(2/). 



(8.; 



By substituting y(x) in (8.8), fti(z) and ft 2 (z) are obtained. Finally it must be tested that 
inequality (6.3) is verified. To enable us to comply with (6.3), the quantity y L may be intro- 
duced as a parameter in 2 (i?). 

Obviously, as the above procedure is purely mathematical, it must be verified "a poste- 
riori'' that the obtained values of U\{x) and n 2 (x) are practically attainable in a given fre- 
quency range. Furthermore, for evaluating n\(x) and n 2 (x), digital computations will be 
generally necessary. Only if the hypothesis n 2 <<ni is satisfied, the use of an analog com- 
puter is practically possible. However, in some special cases, the above procedure can be 
simplified. 

Two different types of (p(y) = <Pi(y)—j<P2(y) functions (to which a favorable E (rj) spectrum 
is related) were chosen, and the corresponding n(x) =rii(x) —jn 2 (x) functions computed by 
means of the above procedure. 

The n(x) functions are traced in figures 6 and 7, under each of which the corresponding 
<p(y) function is listed. The first case was computed by means of an IBM 1620 digital com- 
puter (see sec. 2), and the second one by means of a PACE-TR-10 analog computer (see sec. 2). 

Table 6 gives the power reflection coefficient |p (i?)| 2 percent computed by means of the 
KDE (1.3) and (1.1) plus (1.2), for the two considered cases, and for several values of the 
rate L/\. Under table 6 the numerical values of parameters are, as well, referred. 

In the second case, an initial discontinuity (77 (0) = 1 .25) was considered. 

The results summarized under table 6 seem to be rather poor. But it must be pointed 
out that no interferential effect (at least for the second case) takes place between the two terms 
of (6.12). Better results could probably be obtained by lowering the losses. 




Figure 6. Graphs of the m(x) and n 2 (x) functions 
referring to table 6, 1st column. 
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Figure 7. Graphs of the ni(x) and n 2 (x) functions 
referring to table 6, 2d column. 



Table 6. Power reflection 


coefficient |p (LA)l 2 percent for several values of L/X (absorbing case), and numerical 
values of the parameters 


Functions <p(y) 


1 27T 

In VnM-n!=-4sin y^-y 

, 712 „ , 2ir 

arctg — = B sin -7=7— 2/ 


.4 = 1.4 
B=0.32 

?/l = 6.28 

rp 3.23\0<y<0VL 
6.6iK< v <^ 
= 0.833 


Triangular pulse of height A 
—j rectangular pulse of height 
B 

0<V<Vl 


-.4=0.936 
-6=0.195 

y L =2.u 


Z/X=0.1 


59.9 




76.9 


.2 


31.2 




12.5 


.3 


0.55 




3.3 


.4 


3.82 




0.25 


.5 


1.14 




.02 


.6 


0.65 






.7 


.68 







9. Conclusion 



The possibility of linearizing the basic RDE allows us, in the case of matching', to make 
easy use of both the analysis and synthesis procedures, the last being based on the Fourier 
transform techniques. The extension of the method to the n=n(x) case gives it a large field 
of practical applicability. 

In the case ol absorption, unfortunately, the transform techniques are, generally, applicable 
only to a part of the solution of the RDE, and under rather restrictive conditions : p. must be 
a constant and ni(x) and n 2 {x) must be independent of frequency; owing to the conditions for 
physical realizability, this can be (approximatively) true only in a limited frequency band. 
Only the analysis procedure is, therefore, completely reliable. 

Furthermore, the construction of a stratified layer closely approaching two given rtiix) 
and n 2 (x) functions is a hard task, and. if the number of the strata is not sufficiently high, 
noticeable disagreements with the theoretical values have to be expected. As for the match- 
ing, a "band pass" rather than a "high pass" behavior has to be expected [Franceschetti, 1962]. 
An easy way to design and to construct a broadband absorber may consist of using a single 
material of known and proper electric and magnetic characteristics c(co) and /x(co), and by taper- 
ing it "geometrically/ ' e.g., in the form of dyhedra or pyramids of such dimensions that a 
sufficient number of them is included in a X 2 area. Under the above assumptions, the per- 
formance of this absorber may be considered approximatively equivalent to that of a stratified 
absorber, whose e(x) and p(x) parameters are functions of the "geometrical tapering" and, 
obviously, of the single used material. 
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